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Can robustness against experimental imperfections and noise be embedded into a quantum sim¬ 
ulation? In this paper, we report on a special case in which this is possible. A spin chain can be 
engineered such that, in the absence of imperfections and noise, an unknown quantum state is trans¬ 
ported from one end of the chain to the other, due only to the intrinsic dynamics of the system. We 
show that an encoding into a standard error correcting code (a Calderbank-Shor-Steane code) can 
be embedded into this simulation task such that a modified error correction procedure on read-out 
can recover from sufficiently low rates of noise during transport. 


I. INTRODUCTION 

One of the most promising, and furthest progressed, 
uses of a quantum computer is the quantum simulator, 
wherein the Hamiltonian of one well controlled and un¬ 
derstood quantum system can reproduce the Hamilto¬ 
nian of another system that we wish to study [IH^ ■ The 
challenge for these simulators is the tolerance of noise. In 
principle, this can be done - a universal quantum com¬ 
puter can implement gates consisting of the Trotterised 
Hamiltonian within a fault-tolerant architecture U El- 
However, this digital method of simulation means many 
of the benefits of the original (analog) Hamiltonian sim¬ 
ulation are lost; certainly we can no longer use the com¬ 
paratively easy route of natural Hamiltonian dynamics. 
We seek a middle ground; a method for embedding ro¬ 
bustness into an analog Hamiltonian simulation. Such 
a task appears extremely challenging - even if the sys¬ 
tem’s error channels are as simple as possible (e.g. acting 
locally and independently on each spin), then by the end 
of the simulation of Hamiltonian H, a time t later, the 
error operator O has propagated to The full 

set of errors that we have to adapt to is huge, including 
not only the local errors but also highly non-local ones as 
well. Moreover, without a method for extracting entropy 
from the system throughout the simulation, these tech¬ 
niques can never be scalable [5]; we simply envisage that 
they permit larger, more accurate, simulations than in 
their absence. We report on one special case for which an 
error corrected simulation is possible: perfect quantum 
state transfer. This shows conceptually that additional 
robustness can be imbued upon a quantum simulation, 
and may provide insight that benefits future studies. 

The use of spin chains for transferring a quantum state 
was first proposed by Bose [5] , and refined for perfect ac¬ 
tion in m^- They are intended to reduce the experi¬ 
mental demands of an essential component of the quan¬ 
tum computer - the transport of quantum states between 
distant locations. The concept requires the design, pre¬ 
manufacture, and testing of a device that is made from 
the same technology as the rest of the quantum computer 
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and has a single, fixed function (although ‘routing’ may 
be possible under certain assumptions [HI [T3]). Con¬ 
ceptually, this enables one to expend a lot of effort on 
making the device as accurately as possible to minimise 
errors - since direct control is not required of any of the 
spins in the chain except for the first and last, all the 
rest could (in principle) be isolated from the environ¬ 
ment. Furthermore, it should be comparatively simple 
to experimentally realise this protocol [lU dS] (although 
one has to add remarkably little in order to regain the 
full power of quantum computation m)- 

Inevitably, errors will arise in both the manufacturing 
process, and as noise during the transport [TillTS] . How 
are we to surmount such obstacles? If we allow access to 
a few sites at the beginning and end of a long chain, then 
there is an elegant solution to the task of finding the op¬ 
timal encoding across those spins m in the presence of 
(time independent) identified manufacturing defects. Al¬ 
ternatively, multiple parallel (non-identical) chains may 
be used, at the cost of a heralded, but non-deterministic 
arrival m- Far less is known about dynamically occur¬ 
ring noise, with |19j imposing that errors only occur at a 
restricted number of positions and times, and no true er¬ 
ror correction protocol (for unidentified errors) is known. 

In this paper, we identify the equivalent of ‘local er¬ 
rors’ for a spin chain and investigate error correcting 
codes that can correct for the presence of a small num¬ 
ber of these errors. In the first result of its kind, we show 
that by modifying their error correction procedure, the 
Calderbank-Shor-Steane (CSS) codes HOI El] can be used 
to encode an unknown quantum state into a few sites at 
the beginning of a long chain, and decoded at the op¬ 
posite end, thereby enabling high quality transport of a 
state in the presence of sufficiently low error rates. This 
study is of significant experimental relevance - the free- 
fermion type Hamiltonians which we study have broad 
experimental feasibility , and this work shows how 

they can tolerate (i) perturbative errors in the intended 
coupling strengths (imperfect manufacture), (ii) imper¬ 
fect timing of the state transfer protocol, and (iii) local 
noise which is dominated by one particular type of error. 



2 


II. NOISY TRANSFER CHAINS 

We consider a chain of N qubits with nearest- 
neighbour couplings. The Hamiltonian has the XX form 

N N-1 

i ^ J„(A„A„+1 + Y„Y„+i), 

n—1 n—1 

where Xn,Yn and indicate the standard Pauli X,Y 

and Z matrices respectively, applied to site n (and iden¬ 
tity elsewhere). Such Hamiltonians are particularly ap¬ 
propriate for superconducting qubits [22H21]. The cou¬ 
pling strengths can be selected in many different ways 
to achieve perfect transfer |57j, but the details are irrele¬ 
vant as our constructions are universally applicable - the 
same error correcting code and correction procedure can 
be used on any such chain. Equivalent constructions can 
be made for other one-dimensional free-fermion models, 
such as the transverse Ising model nuns], which is also 
experimentally relevant [3123110], but do not apply to 
the Heisenberg model. For our purposes, it is sufficient 
to know that there exists a time to such that 

for some known phase (j) mi. Here, the A-fold tensor 
product represents the states of consecutive spins on the 
chain; the first is the input spin and the last is the output 
spin of the state transfer process. Eq. Q imposes that for 
any arbitrary initial state Idr/) of the N qubits, after evo¬ 
lution for time to, this state is mirror inverted about the 
centre of the chain, up to the application of controlled- 
phase gates between every pair of qubits [II1I2HI. How¬ 
ever, if a state on a block of spins (such as at either end 
of the chain) has a fixed parity of excitations (number of 
|l)s), then the controlled-phase gates cannot cause that 
block to become entangled with the rest of the system. 
This observation has previously been used to avoid the 
initialisation of any part of the spin chain except where 
the state is input miiMiEni. We seek an encoding for 
an unknown quantum state on the first M spins of a spin 
chain such that, after time to, the input state can be 
recovered from the M spins at the opposite end of the 
chain (the decoding region). In the absence of noise, it is 
sufficient to encode in a state of fixed parity of excitation 
number and use a perfect transfer chain. 

Inspired by the Jordan-Wigner transformation, we in¬ 
troduce the Majorana fermions 

Cjl Z\Z2 . . . Zji—\Xyi Cyj + jV Z\ . . . Zyi—\Yyi, 

whose time evolution can be written as 

2N 

Cn{t) = ^ (m| \n) (2) 

m—1 

where 

h=-Y Z) Hi 


is a 2N X 2N matrix describing the coupling of the 2N 
fermion modes (the tensor product between Pauli Y and 
the matrix Hi is merely a matrix construction and re¬ 
flects no correspondence with a physical division of sub¬ 
sets of qubits), and 

N N-1 

Hi = '^Bn \n) (n| -f ^ J„(|n) (n -|- 1| -I- |n -|- 1) (n|) 

n—1 n—1 

is the Hamiltonian H, restricted to the first excitation 
subspace. The key to this description is that the fermions 
Cn evolve independently of one another (one only has to 
be careful of the ordering of the operators, which con¬ 
tributes the afore-mentioned controlled-phase gates). 

The fermions c„ form a basis that, in principle, any er¬ 
ror (such as Xn) could be described in terms of. However, 
certain errors, such as local phase errors are described by 
only a pair of fermions: = —ic„cjv+„. So, if a phase 

error were to occur on any spin at any time t, then by 
Eq. ([^, at the state transfer time to there would still 
only be two fermions present in the system, and the er¬ 
rors on the output region would consist of no more than 
two operators of the form 

Zn+1-mZn+ 2-M ■ ■ ■ ZN+k-l-M^N+k-M 

Zn+1-mZn+ 2-M ■ ■ ■ ZN+k-l-mYN+k-M 

where k G {1,...M}. These are clearly not the single 
site operators that standard error correcting codes are 
designed to combat. However, observe that whatever er¬ 
ror occurs: (I) There are no more than 2 bit flips {X or 
Y) and, (H) on a site p where there is no bit-flip, there 
is a ^ error only if there are an odd number of bit-flip 
errors on the sites p -I- 1 to A. 


III. CSS CODES 

The CSS codes [13 HII constitute the first known ex¬ 
amples of a quantum error correcting code. An [[n, k, cf|] 
code comprises n physical qubits, encoding k logical 
qubits and is capable of correcting for any [{d — 1)/2J 
single-qubit errors (A, Y or Z). They work by combin¬ 
ing two different classical codes: one for X errors and 
the other for Z errors. After an encoded state has been 
exposed to noise, the original state can be reconstructed 
by first performing a syndrome extraction, in which the 
location of errors is written onto some ancilla qubits, and 
then error correction, in which the ancillas are measured 
and the detected errors are inverted, for each type of er¬ 
ror (A or Z) in turn. 

Assume that we have encoded into a [[M, 1, d]] CSS 
code of distance d > 5, with the additional constraint 
that all the stabilizers and logical operators of the code 
commute with . The perfect mirroring property of 
the perfect transfer system ensures that the code arrives 
perfectly on the decoding region at time to, up to the 
possible controlled-phases that are globally applied. The 
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commutativity of the code with ensures that the 

code space has a fixed parity of excitation number, and 
hence the rest of the chain can be initialised in an ar¬ 
bitrary state and not get entangled with the decoding 
region due to these gates. The internal controlled-phase 
gates are removed by applying controlled-phase gates be¬ 
tween all pairs of qubits fV -|- 1 — M to fV m, returning 
the original code, but updating the errors to 


XN+k-MZN+k+l-MZN+k+2-M ■ ■ ■ Zn 
YN+k-MZN+k+l-MZN+k+2-M ■ ■ ■ Zjsf. 


By performing standard syndrome extraction for the X- 
type errors, we can detect any pair of bit-flip errors. 
Hence, by observation (I) we can detect the location of 
any fermionic operators that lie in the decoding region 
(with the exception of c„CAr+„). This is followed up by 
an error correction step. Obviously, we should apply X 
rotations to the detected error locations as normal. How¬ 
ever, by observation (H) we also apply a Z rotation to 
any spin n that has an odd number of detected bit flips 
on sites 7V-I-1 — M ton — 1. This change in the proce¬ 
dure negates some of the fault-tolerant properties of CSS 
codes |31], but these are irrelevant here. 
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□ 
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Syndrome Extraction: [ ] 
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Logical Phase Errors: [fj] [#>] [ ] 

Syndrome Extraction: □ □ 0 

FIG. 1. (Color online) An error correcting code of 15 qubits 
made up of 3 logical qubits (grouped by shading) in a Z- 
error correcting code, each composed of a 5-qubit A-error 
correcting code. Two bit-flips, and the trailing Z errors, have 
been detected and corrected. Up to two Z errors remain, 
located on the sites of the detected X errors. A single error 
is easily detected and corrected (a). A pair of errors within 
the same logical qubit cancel each other (b). Two errors on 
different logical qubits yield a syndrome that does not match 
the locations of the X errors, flagging the need to correct for 
two Z errors. 


The only errors that we haven’t corrected for are either 
the CnC]sf+n type (a Z on a single site), or the distinction 
between X and Y operators, which were corrected as if 
they were X. Thus, the remaining errors consist of up to 
two Z errors, which can be detected and corrected using 
the Z part of the CSS code in the standard manner. The 
net result is the perfect correction of any single 2-fermion 
error (e.g. Z) that occurs at any time during the evolu¬ 
tion of the system, a far stronger result than the identified 
noise of HU, while by no means contradicting studies of 
information propagation in more noisy scenarios |33j . 

More generally, by selecting an [[M, 1,2fc-|-l]] CSS code 
to encode into (whose stabilizers commute with Z®^), 
we can correct for up to k Major ana fermions arriving 
in the output region. This requires input and output re¬ 
gions of size 0{k^). While such codes are not so well 
studied, [[d^, l,(i]] Shor codes can easily be constructed: 
divide the (P qubits into blocks of d and define logical 
qubits |0)^^ = |-|-)®‘^ and = |—(Hadamard- 

rotated majority vote for Z errors, distance d, where 
|±) = (|0) ± |l))/-\/2), then take those and form them 
into a standard majority vote (i.e. tuned for X errors) 
\Q)l 2 — |0)li and \Yj ^^2 — |1)li ■ This has a fixed parity 
of excitation number provided d is even. However, better 
codes probably exist - in the case of d = 5, a [[19,1, 5]] 
code has been found [33], although this does not have 
the necessary excitation parity condition. Indeed, the 
explicit example that we have simulated in the next sec¬ 
tion suggests that asymmetric quantum error correcting 
codes could become objects of particular interest [35] . 


A. Minimal Working Example 

As presented, the minimal size of the encoding and de¬ 
coding regions is 36 qubits, rendering simulation of the 
full Hilbert space a serious challenge. In order to numer¬ 
ically verify the presented results, we made a number of 
simplifications. Firstly, we limited the length of the chain 
to the size of the encoding region (M = TV), and evolved 
for twice the state transfer time, thereby creating a per¬ 
fect revival of the original state, which also removes the 
requirement for the code to commute with Z®^. Sec¬ 
ondly, we made the observation that although up to two 
Z errors could occur, it is only necessary to use an error 
correcting code that corrects for a single Z error. This 
is because there are two possible cases for what happens: 
either no bit flip error is detected (which means the Z 
error has propagated to a single Z error on some site), or 
two bit flip errors are detected (which means that there 
will be up to two Z errors on those two sites). Obvi¬ 
ously, the only case that a one-Z-error correcting code 
could not implicitly deal with is the instance in which 
there are exactly two Z errors, Fig.j^ However, we know 
that if a one-Z-error correcting code is comprised of 3 
qubits, and there are errors on two of the qubits, the 
syndrome measurement detects an error on the other bit, 
Fig. 0 c). Thus, if we compare this with the syndrome 
information for the X error, we can implement the addi¬ 
tional rule “if a Z error is detected on a logical qubit on 
which no X error was detected, then there were actually 
Z errors on the other two logical bits instead”. As such. 





4 


the Shor-like construction reduces to 15 qubits, render¬ 
ing simulation more feasible. We have implemented this 
m using the standard coupling configuration for perfect 
transfer, J„ = \/n{N — n) and = 0 [ini E]. For 1024 
random samples of a single Z error (chosen randomly to 
occur at any time during the transfer time, and to be 
applied at any random site on the chain), every instance 
was corrected perfectly. 


IV. MORE REALISTIC NOISE 


While a fixed number of errors simplifies the pedagogy, 
of more practical interest is the case of a fixed per-qubit 
error rate. Assume that there is a probability per qubit 
and per unit time of 7 that a two-fermion error occurs 
on the spin chain. For instance, phase errors can be de¬ 
scribed using the Master Equation 


dp 

dt 


N 

-i[H, p] - N-fp + 7 X! 

n—1 


During the time toi one would therefore expect an average 
of 2"fNto fermionic errors to afflict the chain. However, 
we anticipate that only 0{jMtQ) of these are located on 
the decoding region at the moment of error correction. 
To formalise this, consider the intended transformation 
under the perfect state transfer Hamiltonian of c„ 1 —>■ 
CN+i-n and CN+n '—>■ C 2 N+i-n, i-O. we want to keep track 
of errors in the mode for 

n < N. Thus, if we evaluate 

Xn = Tr(pe“*^‘c„e*^*) = Tr(pcAr+i_„), 

where p = is the density matrix in the inter¬ 

action picture, we can evaluate the derivative as 


dXn 

nr = 

N 

\ ^ 7 —iHt ry iHt\ 

T / ^ 'Ll ypc CjV +1 — ) ■ 

m—1 

Using Eq. ([^, we rewrite 

- 2c^ (m| |n) 

-2cN+m {N + m\ |n) 

such that 

^ = - 27 Tr (pe-^^^ ^ (m| |iV + 1 - n) 

\ m=l / 

= -27Xn, 

using Eq. ® again. This leaves a final solution of Xn{t) = 
e“^^‘x„(0L We interpret this as a probability of p = 
i(l — of each fermionic mode having an error. 

If jtQ <C 1, then the error probability is approximately 
P = 7^0 P^r fermionic error. However, these error modes 
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FIG. 2. In the minimal working example, we test faults due 
to timing errors, assessing the probability of successful cor¬ 
rections (by tracking all possible syndrome measurements). 

do not occur independently of one another. Provided 
the expected number of errors ~ 2xMtQ is smaller than 
the number of errors that the code can correct for (~ 
y/M/2), the code can be useful. Thus, we require 7^0 
1/y/M and, furthermore, to scales with at least N if the 
maximum coupling strength of a chain is bounded m- 
We therefore envisage this being applied in finite length 
chains where the parameters can be judiciously chosen to 
be effective, thereby providing regular repetitions of error 
correction if transfer is required over greater distances. 



V. TIMING AND MANUFACTURING ERRORS 

So far, we have shown that a small number of fermionic 
errors can be corrected for, and we have ascribed their 
appearance to noise in a particular basis. However, there 
are two other important mechanisms that can be de¬ 
scribed in this way. The first is a timing error - rather 
than removing the arriving state from the spin chain at 
time to, we accidentally do so at the time to + St [^. In¬ 
stead of the initial state I'l'/) evolving to the target state 
ITt-) = it acquires an error Ex¬ 

panding this for small St (requiring (jtAmax ^ 1, where 
Amax is the largest singular value of Hi) yields an expan¬ 
sion in even powers of the fermionic operators, i.e. larger 
numbers of errors are strongly suppressed such that error 
correction succeeds with high probability. In Fig. we 
show how the minimal working example successfully cor¬ 
rects the majority of cases for a small timing error. Note 
that, unlike previous treatments such as in |38j . here we 
evaluate the probability that the state arrives perfectly, 
not the overlap between the input and output states. 

Similarly, were we to imperfectly manufacture the tar¬ 
get coupling strengths and magnetic fields in H, then 
the perturbation V (which could be time dependent) is 
quadratic in fermions. Now the error in evolution can 
be described by g-dH+v)togiHto and expanded in 
powers of V as 

e-dH+v)to^iHto + r e-^^*Ve^^*dt + 0{tlV^). 

Jo 
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FIG. 3. In the minimal working example, we test faults due 
to coupling errors, altering each uniformly at random in the 
range 7n(l — /) to Jn(l + /), assessing the probability of suc¬ 
cessful corrections and averaging over 1000 different instances 
(black) or the minimum success probability (grey). 

Since the action of preserves the number of 

fermions, this represents an expansion in even powers of 
fermionic operators. Thus, provided Cmaxto ^ 1 where 
Cmax is the largest singular value of the first excitation 
subspace of V, this describes that with high probability 
the number of fermionic errors is small, and these can 
be corrected for via our error correction procedure. The 
efficacy of this procedure is tested for the minimal work¬ 
ing example and depicted in Fig. Of course, a better 
test would be to utilise state transfer over a greater chain 
length to avoid possible confusion with localisation due 
to the errors in the system. 



VI. CONCLUSIONS 

In spite of the fact that CSS codes were designed for 
application in the scenario of local, independent, errors, 
we have shown that they can be retasked to correct for 


the massively correlated errors that typically arise due 
to the intrinsic Hamiltonian dynamics of a spin chain, 
such as those intended to perform the perfect quantum 
state transfer of spin chains. This in turn shows that 
additional robustness against imperfections can be em¬ 
bedded into a fixed-function Hamiltonian evolution such 
as an analog quantum simulation, even if it was depen¬ 
dent on some very specific features of the Hamiltonian 
(the free-fermion structure). 

The major drawback of encoding in this way is that a 
block encoding of size M can tolerate at most 0{VM) 
errors, but a constant per-qubit error rate would re¬ 
quire 0{M) to be tolerated, imposing the requirement 
for rather small error rates. A future direction would be 
to try and improve this through better choice of error cor¬ 
recting code. We are given heart by the surface code [35] 
- this is a CSS code and, although its distance is short, it 
is the case that in the presence of local noise, it is highly 
unlikely that those short error strings that cause failure 
of the code arise. Indeed, the surface code has an error 
correcting threshold consisting of a finite per-qubit error 
rate, exactly as we desire gnjim. However, more work 
would be required - the errors in the present model are 
not independent. Provided these correlations are local, 
it is likely that they could be tolerated, with a mapping 
from the spin chain onto the surface code that preserves 
the locality of the errors. 

The error correction succeeds provided the error op¬ 
erators are well described in terms of a small number 
of fermionic operators. This includes terms such as Z, 
X ® X and V 0 V, as well as Hamiltonian perturba¬ 
tions and timing errors. However, it does not include 
bit-flip errors - a bit-flip on spin n requires 2n — 1 Ma- 
jorana fermions to describe it. It remains an interesting 
question for the future whether either a different error 
correction strategy, or a different class of Hamiltonians, 
permits error correction of all local errors. 
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